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Matchgates are a family of two-qubit gates associated with noninteracting fermions. They are 
classically simulatable if acting only on nearest neighbors, but become universal for quantum com- 
putation if we relax this restriction or use SWAP gates [Jozsa and Miyake, Proc. R. Soc. A 464, 
3089 (2008)]. We generalize this result by proving that any nonmatchgate parity-preserving unitary 
is capable of extending the computational power of matchgates into universal quantum computation. 
We identify the single local invariant of parity-preserving unitaries responsible for this, and discuss 
related results in the context of fermionic systems. 
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I. INTRODUCTION 

There are good reasons to study quantum circuits with 
limited allowed operations. Experimentally, they may 
correspond to circuits implementable with current tech- 
nologies or using a particular physical system. Theoret- 
ically, they can help us understand the minimal require- 
ments for the quantum computational speedup achieved 
in problems such as factoring. When the limited opera- 
tions are not sufficient for universal quantum computa- 
tion, identifying the missing resource can be revealing. 
As an example, the set of single-qubit gates operating on 
a separable initial state can be simulated classically, but 
the addition of any entangling gate to the set makes it 
universal for quantum computation. 

Matchgates are a restricted family of two-qubit gates 
first described by Valiant [1 in the context of graph the- 
ory. Valiant showed that circuits composed only of such 
gates acting on nearest-neighbor qubits were classically 
simulatable. Nearest-neighbor matchgates were shown 
to be equivalent to noninteracting fermions [2], and re- 
main classically simulatable even in the case of adap- 
tive measurements. This is in contrast to noninteracting 
bosons, which allow for universal quantum computation 
with adaptive measurements [3]. 

Many results were obtained by studying the algebraic 
properties of matchgates, with no explicit reference to 
their fermionic nature. For instance, it was shown that 
the computational power of nearest-neighbor matchgates 
is equivalent to that of space-bounded quantum compu- 
tation [4] and of linear threshold gates [5]. It was also 
shown that a universal set can be obtained by relaxing 
slightly the nearest-neighbor condition, allowing match- 
gates to act on next-nearest-neighbor qubits (i.e., one 
qubit apart) or, equivalently, by adding the SWAP gate 
to the set [5]. This algebraic study of matchgates is im- 



portant not only for understanding fundamental proper- 
ties of quantum circuits, but also to study possible ex- 
perimental implementations in different physical systems 
such as linear optics, as proposed in [TJ. 



In this paper we address the question of which re- 
sources are sufficient to extend the computational power 
of matchgates into full quantum computation. We do this 
by studying parity-preserving two-qubit gates, a natural 
generalization of matchgates. We expand on the result 
of [6] by showing that any (nonmatchgate) member of 
this larger set, together with matchgates, is sufficient to 
attain universal quantum computation. Our discussion 
relies on the characterization of two-qubit gates in terms 
of their nonlocal parameters |9] . These are quantities 
invariant by local, single-qubit unitaries; one example is 
the entangling power defined in ; 10j. We show that a 
parity-preserving gate's ability to extend matchgates into 
quantum universality can be attributed to a particular 
nonlocal parameter. This identifies the common char- 
acteristic shared by gates as different as SWAP and CZ 
which enables them to boost the computational power of 
matchgates [51 HI]. 



The paper is structured as follows. In section |II A| we 
define matchgates and parity-preserving gates, recalling 
important known results and demonstrating the univer- 
sality of matchgates + SWAP in a simplified way, and in 
section III Bl we review the characterization of the nonlo- 
cal parameters of two-qubit gates. In section |III A 



we 



suggest an intuitive generalization of the SWAP into a 
continuous family of two-qubit gates, which can provably 
boost the computational power of matchgates. In section 



III B we provide a complete characterization of match- 
gates and parity preserving unitaries in terms of their 
nonlocal parameters, necessary to obtain our main result 
in section III C In that section we show that any non- 
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matchgate parity-preserving gate forms a universal set 
together with matchgates. Sections |IV| and [V] contain a 
discussion of the results, open problems and conclusions. 
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II. REVIEW 

We start by reviewing some basic concepts, definitions 
and known results which will be necessary for our main 
discussion. In section Hi Al we recall the definition of two- 
qubit matchgates and show how they are inserted in the 
larger group of parity-preserving matrices, referring to 
known results regarding simulability and universality of 
these sets. In section [IlB[ we review the characterization 
of the nonlocal parameters of general two-qubit gates. 



A. Parity-preserving unitaries and matchgates 

Let us consider the most general 4x4 unitary ma- 
trix that preserves the parity subspaces ({|00) , |H)} and 
{|10),|01)}) of a two-qubit system. Denoting it by 
G(A, B), where A and B are the 2x2 matrices acting on 
the even and odd subspaces, respectively, we can write it 
as 
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Since G(A,B) ■ G(C,D) = G(AC, BD), the set of all 
unitary matrices of this type forms a subgroup of U(4). 
From this product rule and the unitarity of G(A, B) also 
comes the fact that A and B must be unitary as well. 
With no loss of generality, we can fix the global phase of 
G(A, B) by imposing that det(A)det(_B)= 1. If we also 
impose the nontrivial constraint that det(A)=det(B) — 
± 1, the resulting unitary is known as a matchgate |24j . 
From this point on, parity-preserving (P.P., for short) de- 
notes general unitaries of the form G{A, B), while match- 
gate denotes the special case where the additional deter- 
minant constraint is satisfied. Sometimes, we will also 
refer to "nonmatchgate P.P. unitaries", when det(^4) 7^ 
det(B). 

As mentioned before, the computational power of 
matchgates has been studied in the literature both in the 
context of noninteracting fermions and of the algebraic 
properties of the matrices themselves. The most well- 
known result is due to Valiant lj, who defined match- 
gates in the context of graph theory and proved that 
circuits composed only of nearest-neighbor matchgates 
can be classically simulated. This result was reproduced 
in the context of fermions by Terhal and DiVincenzo 
[2J, who showed that circuits composed only of nearest- 
neighbor matchgates were mathematically equivalent to 
systems of noninteracting fermions. They also showed 
that the simulability result holds even in the presence 
of adaptive measurements, in contrast to noninteracting 
bosons, for which adaptive measurements allow for uni- 
versal quantum computation [3]. Subsequently, several 
results were obtained regarding how to increment this 
set of operations to attain universal quantum computa- 
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FIG. 1: Universal set of gates in terms of matchgates and 
SWAP: (a) Arbitrary single-qubit gate A and (b) CZ gate. 

tion, for example by adding certain simple interactions 
[11] or allowing for charge measurements (12| . 

In this work, we are interested in the algebraic proper- 
ties of matchgates. In particular, we address the question 
of which unitaries can single-handedly boost the compu- 
tational power of (nearest-neighbor) matchgates to that 
of full universal quantum computation. One such exam- 
ple is the SWAP gate, since it replaces the need for long- 
range interactions, and it is known [2J [T3] that match- 
gates acting on arbitrary pairs of qubits form a universal 
set. Even more, it has been shown in [BJ that one does 
not need to fully remove the nearest-neighbor constraint, 
and only relaxing it slightly to allow for next-nearest- 
neighbor interactions (i.e., between qubits one line apart) 
is sufficient to achieve (encoded) universality. 

Since our work is heavily based on the construction 
proposed in [BJ, we now reproduce a simplified demon- 
stration of this result. In [BJ, each logical qubit is encoded 
in four physical qubits. Here, we will use an encoding 
which is more economical (each logical qubit is encoded 
in two physical qubits), and which is equivalent up to 
some technical reasons which need not concern us. We 
encode each logical qubit in the even parity subspace of 
two physical qubits: 

|0) L = |00), 

ll)i = lll>. 

Since matchgates act separately on parity subspaces, 
any single-qubit gate A is implemented in the encoded 
qubit simply by the matchgate G(A, A), as shown in Fig. 

[Qi. 

Besides arbitrary single-qubit gates, for universality it 
is sufficient to implement one of the standard maximally 
entangling gates. One choice specially suitable for this 
construction is the controlled- Z [CZ = diag(l, 1, 1, — 1)], 
depicted in Fig. [lp. Let us label the physical qubits that 
encode the first logical qubit with numbers 1 and 2, and 
similarly with the second logical qubit (physical qubits 
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3 and 4). In view of the encoding, it is simple to see 
that a CZ between the logical qubits can be implemented 
by a single CZ between physical qubits 2 and 3 only. 
This is the simplest choice which preserves the encoding 
while still avoiding the need for long-range interactions 
(compare with the CNOT gate, for example, which would 
have to act on qubits 3 and 4 controlled by the state of 
either qubit 1 or 2). 

It can be easily checked that the circuit in Fig. [T|b, 
implementing 

CZ 23 = G(H, H)G(X, X) SWAP G(H, H), (3) 

has the desired effect and thus completes the universal 
set (in the encoded space) , having used only matchgates 
and the SWAP gate. Note that qubit initialization and 
measurements are not an issue here, as both preparation 
and measurement in the logical computational basis are 
trivially implemented by preparation and measurement 
in the physical computational basis, and there is no need 
for previously entangled or otherwise inaccessible states. 
In short, this gate-set can efficiently simulate any quan- 
tum computation. 

There are other examples in quantum information of 
a single gate uplifting a classically simulatable set into 
universal quantum computation. For instance, circuits 
with nonentangling gates only are classically simulat- 
able, while adding a single entangling gate allows for full 
quantum computation. Another example is the Toffoli 
gate, which by itself results in classically simulatable cir- 
cuits (in fact, Toffoli is universal for classical computa- 
tion), but adding the Hadamard gate promotes the set 
to quantum universality [H]. However, the matchgate 
+ SWAP result is also surprising as the added resource 
seems to be a very trivial one (the simple swapping of 
qubits), as opposed to the resources in the other cases, 
which are considered to be fundamentally quantum in na- 
ture (entanglement and quantum superposition, respec- 
tively). Also of notice is that, in contrast with the Toffoli 
gate, matchgates are expected to be much weaker even 
than full classical computation, as it can be shown that 
any boolean function calculated by a circuit composed of 
nearest-neighbor matchgates is trivial, in the sense that it 
depends only on one bit of the input [5]. In this sense, the 
SWAP gate is bridging the gap between a sub-classical 
complexity class and full quantum computing. 

In [B] , the role of the SWAP gate was seen as to provide 
the long-range interactions missing in nearest-neighbor 
matchgates. In this case, we could expect other gates 
that exchange the role of the qubits in the computational 
basis to do just as well, such as iSWAP [= G(I, iX)], and 
fermionic SWAP [= G(Z,X)]. It is easy to check, how- 
ever, that these two gates are matchgates, and thus can- 
not confer universality to matchgates as SWAP can. An- 
other aspect in which they differ from the SWAP is that 
they are entangling. The fermionic SWAP, for example, 
is a perfect entangler (it can be written as CZ-SWAP). 
Their entangling action might create unwanted correla- 
tions that corrupt the computation, and be the reason 



why neither gate can play the role of SWAP in this con- 
text. 

We will clarify this question by pinpointing precisely 
the characteristic that a parity-preserving gate must have 
in order to form a universal set together with matchgates. 
Before we address this question, we must turn to the 
characterization of two-qubit unitaries in terms of their 
nonlocal parameters. 



B. Nonlocal parameters 

Let us begin this section by stating a result which will 
be fundamental for the discussion that follows: 

Theorem II.l JH \Wj Any two-qubit gate U G SU(4) 
can be written as 

U ={U 1 ®U 2 )U NL {V 1 ®V 2 ) 

= (Ui ® u 2 y(*x®x + bY®Y+cz®z) (Vi ^ y 2) (4) 

where Ui and V% are single- qubit gates on qubit i, and a, 
b and c are real parameters in the interval [0, 5). 



In other words, Theorem II.l states that, out of the 



15 parameters that define a general 4x4 unitary matrix 
(up to a global phase), only 3 are truly nonlocal (in the 
sense that they require some interaction between the two 
qubits), while the remaining 12 can be implemented by 
two sets of single-qubit gates, to the left and to the right 
of this nonlocal "core" . 

The set of nonlocal parameters described above is con- 
venient because it gives an explicit way to parametrize 
two-qubit unitaries, but it does not seem to offer much 
physical insight. One reason is the non-uniqueness of the 
decomposition, as we can implement permutations of the 
type a O c using only single-qubit gates (in this case, 
H® 2 ). This generates an ambiguity if one tries to char- 
acterize unitaries in terms of the triplet {a, b, c}, since 
any unitary gate U can be associated with a (discrete) 
number of different sets of parameters, related to each 
other by permutations. One usual way of dealing with 
this is choosing a > b > c or a similar condition, which 
eliminates the ambiguity by fixing one among all possible 
sets for a particular matrix (alternatively, for a geomet- 
rical approach in terms of what is known as the Weyl 
chamber, see [§])• As we will see in the next sections, 
however, for the sets of matrices we are interested in this 
is not an issue, and we will make no such restriction. 

While {a, b, c} are not, strictly speaking, local invari- 
ants, one can always write true local invariants as (sym- 
metric) functions of these parameters. One such example 
is the set of invariants defined by Makhlin in [5] , as shown 
in [5]. Another example, in which we will be particularly 
interested in the discussion that follows, is the entangling 
power defined by Zanardi et al. in |10j . 

The entangling power e p (U) of an unitary U is defined 
as the average entanglement generated by the action of 
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U on all product states IV'i) <8> IV^) [TOj : 



the set of nearest-neighbor matchgates in order to create 
a universal set. In section III A we define a continu- 



(5) 



where the bar denotes average with respect to some 
probability distribution p(ipi, "02 )• Simple manipulations 
show that this is equal to 



s p (U) = 2 tr [U® 2 Q P U^ 2 P U ] , 



(6) 



where P± 3 is the projector on the antisymmetric subspace 
of qubits 1 and 3 (notice that the trace is taken over a 
doubled Hilbert space), while 



n p = J dn{i>i,,h){\<h) m ® (H \Tp 2 )) m , 

and d/j, denotes the measure over the space of product 
states induced by probability distribution p(ipi,ip2). 

It can be shown that, if the average is taken over the 
uniform distribution, the entangling power is both lo- 
cal invariant and SWAP invariant (that is, it remains 
the same if U is multiplied by the SWAP or by single- 
qubit gates on either side), which, as previously men- 
tioned, suggests that it must be written as some function 
of {a,b,c}. If fact, simply by plugging Eq. Q in the 
above definition (and rescaling the value of e p (U) so it 
goes from 0, for local gates and SWAP, up to 1, for per- 
fect entanglers) we find that 

e p (U) = 1 - cos 2 2a cos 2 26 cos 2 2c - sin 2 2a sin 2 26 sin 2 2c. 

(7) 

It is clear from the expression above that any local 
gate has e p (U\ <S> U2) = 0. It can also be seen that 
e p {CNOT) = 1, since for the CNOT {a,b,c} = {f ,0,0}, 
and that e p (SWAP) = 0, since for the SWAP we have 
{£ , j , j }, all results which were expected. We also see 
that the above expression is invariant under permutations 
of {a, b, c}, another previously anticipated feature. 

Recall that the set of single-qubit gates together with 
any entangling gate is universal for quantum computa- 
tion. We can check that a gate is entangling by calcu- 
lating the non-local invariant defined above. Can a sim- 
ilar characterization be given in the context of match- 
gates? The SWAP gate, which seems to be responsible 
for boosting their computational power, cannot be imple- 
mented only by single-qubit gates, and thus must clearly 
have some nonlocal property (even though it is not en- 
tangling). We will show in the next section that this 
property can be characterized in terms of the {a, b, c} 
parameters above, and that it is fundamental for grant- 
ing computational power to matchgates. 



III. EXTENDING MATCHGATES WITH 
PARITY-PRESERVING UNITARIES 

In this section, we present our main result, regarding 
which parity-preserving (P.P.) unitaries can be added to 



ous family of P.P. gates that interpolates between the 
SWAP and iSWAP, in order to understand why SWAP 
can boost nearest-neighbor matchgates to quantum uni- 
versality whereas iSWAP cannot. In section |III B| wc 
characterize an arbitrary P.P. gate in terms of its nonlo- 
cal parameters, identifying both the subgroup of match- 



gates and the continuous family of section III A in this 
context. Finally, in section [ill C| we prove the main the- 
orem by determining under which conditions a specific 
P.P. gate can replace the SWAP in the universal set of 



section II A and explicitly constructing a universal set 
whenever these conditions arc satisfied. 



A. Interpolating between SWAP and iSWAP 

Let us begin by addressing a question raised at the end 
of section II A| why is it that, in the context of nearest- 
neighbor matchgates, the SWAP gate is sufficient to gen- 
erate a universal set, but other gates similar to it (such as 
iSWAP or fermionic SWAP) are not? In order to answer 
this question, let us replace the SWAP by a continuous 
family that interpolates between it and iSWAP, given by 



G(I, e lT X) 
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(8) 



where r goes from (SWAP) to tt/2 (iSWAP). This is 
clearly a P.P. gate [cf. Eq. Q], but it is also a matchgate 
whenever t — tt/2 (i.e. only in the case of the iSWAP). 
Thus, we expect the gate that arises from this substi- 
tution in the circuit of Fig. [T]d to complete the univer- 
sal set (together with the single-qubit gates from Fig. 
[Ija,) for one extreme value of r, and to fail to do so for 
the other extreme. What about intermediate values? In 
those cases, the resulting gate implemented in the logical 
qubits is given by G{H, H)G(X, X)G(I, e iT X)G(H, H) 
or, explicitly, by: 



/ 1 

e iT 

e lT 

Vo 0-1 



(9) 



This matrix is diagonal, and thus preserves the encod- 
ing [cf. Eq. Q], which is a requirement for it to be used 
in an encoded universal set. Simple calculations show us 
that the three nonlocal parameters [cf. Eq. Q] for this 
gate are {0,0, | — |} and so, by Eq. Q, we find that 
its entangling power is cos 2 r. This confirms our previ- 
ous intuition: the entangling power is maximum when 
we use the SWAP, but goes to as the gate used ap- 
proaches the iSWAP, in which case any circuit composed 
of these gates must become classically simulatable. The 
simulability of the resulting circuit can be understood in 
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two ways: either by remembering that, in this case, the 
whole circuit is composed of matchgates, or by noticing 
that, when the two-qubit gate creates no entanglement, 
any circuit acting on the logical qubits can be expressed 
in terms only of single-qubit gates. 

However, perfect entanglers are not necessary to form 
a universal set |17j . If arbitrary single-qubit gates are 
available, a gate that creates any (nonzero) amount of 
entanglement is sufficient for universal quantum compu- 
tation, with at most a polynomial overhead with respect 
to the set with perfect entanglers. In view of this, we see 
that the G(I, e lT X) gate above is sufficient to achieve uni- 
versal quantum computation together with matchgates if 
and only if t ^= ir/2 - or, equivalently, if and only if it is 
not a matchgate. 

This shows that there is a continuous family of gates 
akin to the SWAP which can be used to create a universal 
set. Furthermore, the inability of iSWAP or fermionic 
SWAP to provide the extra computational power cannot 
be attributed to the fact that they are also entanglers, 
since the family G{I, e 1T X) includes members with all 
possible values of entangling power. While it is true that 
the most efficient gates in this set are the least entangling 
(which is evident by the particular cases of the SWAP 
and iSWAP), it will become clear in the next sections 
that this is just an attribute of this particular family, 



and there is no such relation for P.P. gates in general. 

This shows that the parameter r in the parametriza- 
tion above seems fundamental for the gate's ability to 
generate a universal set together with matchgates. In 
what follows, we shall relate r to the nonlocal parameters 
defined in the previous section and to the relative phase 
between the submatrices of a P.P. gate (section II A I, in- 
vestigating its importance to the computational power of 
general P.P. matrices when associated with matchgates. 



B. Nonlocal parameters of P.P. matrices 

Let us start by obtaining a convenient parametrization 
of general P.P. matrices. To this end, let us recall that 
any single-qubit unitary matrix can be parametrized as 



cos 9e i( - p+a ^ i sin 9e l ^ + ^ 



where 9, (3, a and 7 are real parameters in the interval 
[0, 2w] . Notice that the determinant of this matrix, e 2l/3 , 
depends only on /3, which will be crucial for the following 
discussion. By parameterizing A and B this way, the 
most general P.P. matrix can be written as 
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(10) 



in terms of 7 real parameters: {9, a, 7, cj>, fi, v, /?}. Notice that, while A and B are each described by 4 parameters, 
G(A, B) is defined only up to global phase, which means we can choose the determinant parameter suitably such that 
det(A) = det(B)- 1 =e 2 ^. 

How much nonlocal freedom remains in the 7-parameter family above? In other words, how does the P.P. restriction 
constrain the nonlocal parameters of G(A, B)l To answer this, let us explicitly write the nonlocal "core" defined in 
Eq. @: 



U, 



NL 



( cos(a - b) e lc i sin (a - b) e ic \ 

cos(a + b) e~ lc i sin (a + b) e~ lc 

i sin (a + b) e~ lc cos (a + b) e~ lc 



y i sin (a — 6) 











(11) 



cos(a — b)e lc J 



Comparing Eq. (10) and Eq. (11), we see that Unl lent to fixing one of these nonlocal parameters (c = 0). 



is, itself, a P.P. gate. Not only that, but simple inspec- 
tion shows us that the nonlocal parameters of G(A, B) 
are given by {—^-, ^y^,/3}, which points us to two im- 
portant facts. First, that general P.P. gates can have any 
value of the nonlocal parameters (and, consequently, of 
any local invariant derived from them) , which means that 
any matrix in U(4) is locally equivalent to a P.P. gate. 
Second, that the matchgate condition (/3 = 0) is equiva- 



An important remark must be made here, regard- 
ing the previously mentioned ambiguity of the triplet 
{a, b, c}. In the previous section, we mentioned that 
permutations between these parameters could be imple- 
mented using only local gates, which meant these param- 
eters were not local invariants per se (and any true in- 
variant, such as the entangling power in Eq. (J7|, must be 
symmetric under these permutations). Notice, however, 
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that the local gates which implement these permutations 
are, in general, not parity-preserving (the Hadamard 
gate, for example). In fact, as will become clear below, 
all single-qubit P.P. gates commute with the Z ®Z gate. 
This means that the P.P. condition distinguishes param- 
eter c from the other two in a qualitative way, by making 
it a true local invariant (and so our characterization of 
matchgates as the particular case of P.P. gates with c = 
is meaningful). 

Now compare Eq. (10) and Eq. (11) again. Three 



of the independent parameters of G{A, B) have been ac- 
counted for, in terms of {a,b,c}. The remaining four 
parameters are given by individual phases, and can be 
achieved by multiplying whole rows and columns by 
phases (i.e., by applying Z- rotations to the left and to 
the right). In fact, a simple calculation shows that 

G(A,B) = [R z { X i)®Rz{x2)]Unl [Rz{&)®Rz{&)], 

(12) 

with an appropriate choice of xi> X2, £i an d £2 in terms 
of a, 7, fi and v. Now, notice that these Z-rotations can 
be rewritten as 

R z (xi) ® Rz(X2) = diag ( e ^(xi+X2) ; e ,( Xl -X2) j 

gi(X2-Xl) g*(-Xl-X2) 

= G(R z (xi+X2),R*(xi-X2)), 



which is a matchgate. Together with Eq. (12), this leads 



to a simple result which will be important later on: 
Lemma III.l Any P.P. gate can be written as 

G(A, B) = G(R z (Ti),R z (T 2 ))e l ( aX ® x+bY ® Y+cZ ®^ 
x G(R z (n),R z (n)). 



(12) 



The importance of this result is that it offers a de- 
composition for P.P. gates in terms of matchgates and a 
single nonmatchgate P.P. gate. Since we are interested 
in the most general P.P. gate with the ability to enhance 
the computational power of matchgates, we will be able 
to disregard the Z-rotations in the above decomposition. 
Being matchgates themselves, they cannot offer any com- 
putational power beyond that of the whole matchgate 
group, and it will be sufficient to consider only the non- 
local term in the middle [i.e. Unl in Eq.Q]. 



C. Main result 

We are now in condition to generalize the result of sec- 
tion . 



Ill A by searching for the most general G(A, B) gate 



to replace the SWAP alongside matchgates in the univer- 
sal set. First, let us consider how far we can generalize 
the circuit in Fig. [TJj. 

In view of the encoding given in Eq. ^ we see that, in 
order to consider gates acting only on the pair (2, 3), the 



resulting gate must be diagonal. Any nondiagonal com- 
ponent would change the parity of one of the pairs (1,2) 
or (3,4), thus disrupting the encoding. By the property 
that G(A, B) is block diagonal, we see that, if both A 
and B satisfy 

H X AH = e ic R z {a x ), 
H X B H = e^ tc R z (a 2 ), 

then the resulting gate is 

G(e ic R z (cr 1 ), e~ ic R z (a 2 )), 

which is the most general diagonal gate. 

This means that both A and B must be, at most, X- 
rotations (multiplied by a phase), and it is clear that 
the family considered previously (which interpolates be- 
tween SWAP and iSWAP) is a particular case of this. 
The most general case, thus, is a matrix of the form 
G(e lc R x (9),e~ ic R x (</>)), given explicitly by 



/ cos(0) e ic 



V % sin (9) e l 





cos(4>) e~ 
i sin (0) e" 






i sin (0) e" 
cos((p) e~ 




i sin (9) e l 



cos(0) e lc 



\ 



= / 
(13) 

which has the exact form of the nonlocal gate Unl, if we 
take a — b = 9 and a + b = <p. 

Thus, we conclude that the most general gate which 
can replace the SWAP in the circuit in Fig. [T]d while 
still preserving the encoding is of the form Unl — 
exp [i(aX ® X + bY (g) Y + cZ (g) Z)\ which, by the re- 
mark at the end of the previous section, also happens 
to be the most general P.P. gate we need to consider in 
order to characterize the ability of any P.P. gate to boost 
the computational power of matchgates. 

By replacing the SWAP with the gate in Eq. 
ph above, the resulting gate, G(H, H)G(X, X)U N l 



x G(H, H), is given explicitly by (in terms of {a, b, c}): 



/ gi(a-6+e) 





V 





a i(a+b- 









i{— a — 6— c) 









— a+b+c) 



\ 



)) 



(14) 

Again, in the same fashion as in section [ill A[ we must 
determine the gate's entangling power. As long as the en- 
tangling power is nonzero, this gate completes the univer- 
sal set together with the single-qubit gates of Fig. [IJl By 
obtaining the nonlocal parameters of the above gate and 
using expression (JT]), we find that its entangling power is 
sin 2 (2c). Again, as before, this means that the gate cre- 
ates some entanglement as long as c 7^ - that is, as long 
as it is not a matchgate. We now have all ingredients for 
the proof of our main result: 

Theorem III.l Let G(A,B) be any parity-preserving 
unitary, and M be the set of all matchgates + G(A,B). 
Then, if G(A, B) is not a matchgate [i.e. if det(A) 7^ 
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T 
I 



G(H,H) 



R"z«i) 



R Z (U 



G(A,B) 



R z (x 1 : 



R' z (xJ 



G(X,X) G(H,H) 



M of matchgates + G(A, B) is universal for quantum 
computation or it generates only classically simulatable 
circuits. This did not need to be so, as there are re- 
cent results about families of quantum circuits which are 
expected to be neither universal for quantum computa- 
tion nor classically simulatable. Examples include the so- 
called IQP circuits [TS] and linear optics without adaptive 
measurements I19I. 



FIG. 2: Circuit where the SWAP is replaced by G(A,B) and 
some additional matchgates. 



IV. DISCUSSION 



det(B)], a universal set can be created from gates in M 
acting only on nearest-neighbor qubits. 

Proof: Let us give a constructive proof by supposing one 
has the ability to use nearest-neighbor matchgates when- 
ever necessary, in addition to a particular G(A, B) gate. 



1. Encode logical qubits as |0) L = 1 00) and |1) 



111) 



2. Obtain the parameters for matrix G(A, B) such that 
it is written in the form 

RAxi) ® R z (x2)e l ^ x+bY ^ Y+cmz) R z (^) ® 

3. Use the circuit in Fig. [T^, to implement any single- 
qubit gate; 

4. Modify the circuit of Fig. [Tp such that the Z rota- 
tions cancel out, and all that is left is 

G(H, H)G(X, xy (aX ® x+bY ® Y+cZ ® z) G{H, H); 

this is an allowed procedure, since the gates required 
to modify the circuit are matchgates themselves. The 
final circuit can be seen in Fig. [2] 

5. The modified circuit now implements the unitary 



/ gi(a-6+c) 



V 



o 

i{a+b—c) 









_gi(— a— 6-c) 




_ e i(-a+b+c)^ J 



which generates entanglement sin 2 (2c). By a known 
result |17j , we know that any entangling matrix, along 
with arbitrary single-qubit gates, is universal for quan- 
tum computation. So, the gate implemented by the 
circuit in Fig. [2] completes the universal set as long as 
c ^ or, in other words, as long as it is not a match- 
gate. 



Note that this result can be complemented by the 
well-known classical simulability of matchgate circuits 
[2]. This means that, at least when restricted to parity- 
preserving nearest-neighbor operations, the jump from 
(sub-)classical computation to full universal quantum 
computation is abrupt, in the sense that either the set 



So far, we have restricted our attention to P.P. gates, as 
they are a natural generalization of matchgates. Which 
other SU(4) unitaries can, together with matchgates, 
form a universal set? It is easy to see that matchgates 
+ arbitrary single-qubit gates form a universal set, since 
matchgates include many perfect entanglers [2J. We also 
know that the Hadamard gate alone is sufficient for this 



purpose since, as mentioned in section II B conjugating a 
P.P. unitary by H® 2 implements the permutation a f-> c. 
Thus we could take a matchgate with a / and, us- 
ing the H gate, construct a gate with c ^ satisfying 
the conditions of Theorem Clearly, other ii-like 

gates which implement similar permutations (i.e. a <R- c 
and b -f-> c, but not a <-> b) can be used in the same 
way. It can also be shown easily that such gates are the 
only single-qubit gates which take matchgates into gen- 
eral P.P. gates by conjugation. While the ability of H to 
extend the power of matchgates appears very naturally 
in the present construction, it is not particularly sur- 
prising: matchgates contain all single-qubit ^-rotations 
which, when added to the H gate, can generate arbitrary 
single-qubit gates, completing a universal set even with- 
out encoding. 

Our results used minor modifications of the encoding 
proposed in [5] and rely heavily on it. Since we encoded 
each logical qubit in the parity of the physical qubits, we 
considered only P.P. gates or, alternatively, operations 
which are not P.P. themselves, but which preserve a cer- 
tain subset of P.P. unitaries (i.e. conjugation by H® 2 ). It 
is an open question whether different encodings may help 
identify different unitaries which generate a universal set 
together with matchgates. 

Another question we may ask is whether there are 
other interesting subsets of P.P. gates which can be shown 
to be either classically simulatable or quantum univer- 
sal. Since matchgates are P.P. gates with c = 0, two 
other natural subsets to consider would be P.P. gates 
with a = or b = 0. Let us consider this latter case, 
and denote by Mb the subset of all P.P. gates with 6 = 0. 
Notice that, in Figs, [iji and [2] only G(A, A) matrices 
are used, along with an arbitrary G(A, B) and local Z- 
rotations. Inspecting Eq. ( [To| ) we see that, in addition 
to having c = 0, G(A, A) has b = (i.e. = 6). This 
means that the present construction already uses mostly 
gates from Mb- By choosing G(A,B) as any gate from 
Mb with c ^ 0, we satisfy the conditions of Theorem 



III. 1 and thus complete a universal set. This proves 
that P.P. gates with b = 0, or equivalently with a = 0, 
are already universal for quantum computation, in sharp 
contrast with the c = case. Thus, the theorem ac- 
tually ascribes a special status for parameter c, which 
we already anticipated in section [ill B| by showing it was 
qualitatively different from the others. 

The above result can be understood easily in the con- 
text of the fermionic nature of matchgates. As can be 
seen in [2j 16] , matchgates are the group generated by the 
set {X <g> X, Y ® Y, 1 <g) Z, Z <8> 1, X <g> Y, Y <g> X}. 
Using the Jordan- Wigner transformation, these matrices 
are shown to be equivalent to Hamiltonians quadratic in 
creation and annihilation operators, describing the evo- 
lution of noninteracting fermions. On the other hand, 
the Z ® Z matrix, related to parameter c in our con- 
struction, is equivalent to a Hamiltonian quartic in cre- 
ation/annihilation operators, and so it describes an in- 
teraction. Thus, the break in the permutation symmetry 
imposed by the P.P. condition, which is not obvious in 
the algebraic approach, turns out to have a very physical 
meaning in the fermionic approach. 

In fact, in [TT], Bravyi and Kitaev have shown that 
the exp{ijZ ® Z} interaction can be used to implement 
a CZ, which is then used to generate a universal set 
with fermionic modes. Our result is similar in spirit to 
this one, but also more general as we provide an explicit 
construction where any nonmatchgate G(A, B) matrix 
is sufficient to generate a universal set. This alterna- 
tive approach can be useful if one is attempting to im- 
plement circuits of P.P. gates in other physical systems, 
rather than fermions. In this case, it is not obvious what 
property gates such as CZ and SWAP have in common, 
except for their ability to provide matchgates with full 
quantum computational power. The power of the SWAP 
is attributed to the long-range interaction it simulates, 
whereas the CZ is used to generate entanglement bet ween 
the encoded qubits (see Fig. [T]d). Thus, Theorem III. 1 



connects both results by identifying the property these 
gates share, showing that it is precisely the resource miss- 
ing in the group of matchgates and providing an explicit 
construction for any other gate of this kind. 

Other interesting subsets of P.P. gates for which simi- 
lar results are known are the unitary evolutions generated 
by Heisenberg isotropic (Hi = X(g> X + Y&Y + Z&Z) or 
anisotropic (H a = X ® X + Y ® Y) interactions [151120] . 
These sets are of great physical interest, as they provide 
examples where a quantum computer could be built out 
of a single interaction Hamiltonian. Explicit construc- 
tions for said universal sets can be found, respectively, in 
[2T] and [22] . One fundamental difference between these 
results is that the isotropic interaction considered in [21] 
acts only on nearest-neighbors, while the anisotropic in- 



teraction in [22] acts on next-nearest-neighbors. Further- 
more, it has been shown that the anisotropic interaction 
is not universal while acting only on nearest-neighbors, 
meaning that distant interactions are actually necessary 
in this case [23]. These two results are clearly consistent 
with our own, as any unitary evolution generated by H a 
must be a matchgate, while Hi generates nonmatchgate 
P.P. unitaries. It is important to point out, however, that 
neither result implies or is implied by our own. 



V. CONCLUSION 

In [6] it was shown that the SWAP gate can be used 
to extend the computational power of nearest-neighbor 
matchgates to full quantum computation. We have ex- 
panded on this result by showing that the SWAP can be 
replaced by any parity-preserving (P.P.) two-qubit gate 
which is not a matchgate itself. We have also character- 
ized arbitrary parity-preserving gates in terms of their 
three nonlocal parameters {a, b, c}, showing that match- 
gates form a particular subgroup of these where c = 0. 
Thus, parameter c is seen to be qualitatively different 
from the other two. This is supported by the fact that 
the set of P. P. gates with b = is universal for quantum 
computation (which is true also for gates with a = 0). 

We have also discussed the relationship between 
our work and other known results in the context of 
fermions. While matchgates are equivalent to noninter- 
acting fermions, the parameter c is known to be associ- 
ated with certain interactions, providing an underlying 
physical interpretation to its special status. Our result, 
however, is more general in the sense that it does not 
rely on the fermionic nature of matchgates, being also 
suitable for implementations in other physical systems. 

We hope that the characterization of unitaries in terms 
of their nonlocal parameters will prove fruitful in under- 
standing further aspects of the separation between classi- 
cal and quantum computation. On a more practical note, 
we believe it could open the way to the construction of 
different universal sets of quantum operations, perhaps 
more adequate for particular experimental quantum com- 
puter implementations. 
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